Multiomics modeling of the immunome, transcriptome, microbiome, proteome and metabolome adaptations during human pregnancy.

MotivationMultiple biological clocks govern a healthy pregnancy. These biological mechanisms produce immunologic, metabolomic, proteomic, genomic and microbiomic adaptations during the course of pregnancy. Modeling the chronology of these adaptations during full-term pregnancy provides the frameworks for future studies examining deviations implicated in pregnancy-related pathologies including preterm birth and preeclampsia.ResultsWe performed a multiomics analysis of 51 samples from 17 pregnant women, delivering at term. The datasets included measurements from the immunome, transcriptome, microbiome, proteome and metabolome of samples obtained simultaneously from the same patients. Multivariate predictive modeling using the Elastic Net (EN) algorithm was used to measure the ability of each dataset to predict gestational age. Using stacked generalization, these datasets were combined into a single model. This model not only significantly increased predictive power by combining all datasets, but also revealed novel interactions between different biological modalities. Future work includes expansion of the cohort to preterm-enriched populations and in vivo analysis of immune-modulating interventions based on the mechanisms identified.Availability and implementationDatasets and scripts for reproduction of results are available through: https://nalab.stanford.edu/multiomics-pregnancy/.Supplementary informationSupplementary data are available at Bioinformatics online

Expression of specific inflammasome gene modules stratifies older individuals into two extreme clinical and immunological states

Low-grade, chronic inflammation has been associated with many diseases of aging, but the mechanisms responsible for producing this inflammation remain unclear. Inflammasomes can drive chronic inflammation in the context of an infectious disease or cellular stress, and they trigger the maturation of interleukin-1β (IL-1β). Here we find that the expression of specific inflammasome gene modules stratifies older individuals into two extremes: those with constitutive expression of IL-1β, nucleotide metabolism dysfunction, elevated oxidative stress, high rates of hypertension and arterial stiffness; and those without constitutive expression of IL-1β, who lack these characteristics. Adenine and N4-acetylcytidine, nucleotide-derived metabolites that are detectable in the blood of the former group, prime and activate the NLRC4 inflammasome, induce the production of IL-1β, activate platelets and neutrophils and elevate blood pressure in mice. In individuals over 85 years of age, the elevated expression of inflammasome gene modules was associated with all-cause mortality. Thus, targeting inflammasome components may ameliorate chronic inflammation and various other age-associated conditions

Systemic Immunologic Consequences of Chronic Periodontitis

Chronic periodontitis (ChP) is a prevalent inflammatory disease affecting 46% of the US population. ChP produces a profound local inflammatory response to dysbiotic oral microbiota that leads to destruction of alveolar bone and tooth loss. ChP is also associated with systemic illnesses, including cardiovascular diseases, malignancies, and adverse pregnancy outcomes. However, the mechanisms underlying these adverse health outcomes are poorly understood. In this prospective cohort study, we used a highly multiplex mass cytometry immunoassay to perform an in-depth analysis of the systemic consequences of ChP in patients before (n = 28) and after (n = 16) periodontal treatment. A high-dimensional analysis of intracellular signaling networks revealed immune system–wide dysfunctions differentiating patients with ChP from healthy controls. Notably, we observed exaggerated proinflammatory responses to Porphyromonas gingivalis–derived lipopolysaccharide in circulating neutrophils and monocytes from patients with ChP. Simultaneously, natural killer cell responses to inflammatory cytokines were attenuated. Importantly, the immune alterations associated with ChP were no longer detectable 3 wk after periodontal treatment. Our findings demarcate systemic and cell-specific immune dysfunctions in patients with ChP, which can be temporarily reversed by the local treatment of ChP. Future studies in larger cohorts are needed to test the boundaries of generalizability of our results

Single-cell Analysis of the Neonatal Immune System Across the Gestational Age Continuum

Although most causes of death and morbidity in premature infants are related to immune maladaptation, the premature immune system remains poorly understood. We provide a comprehensive single-cell depiction of the neonatal immune system at birth across the spectrum of viable gestational age (GA), ranging from 25 weeks to term. A mass cytometry immunoassay interrogated all major immune cell subsets, including signaling activity and responsiveness to stimulation. An elastic net model described the relationship between GA and immunome (R=0.85, p=8.75e-14), and unsupervised clustering highlighted previously unrecognized GA-dependent immune dynamics, including decreasing basal MAP-kinase/NFkB signaling in antigen presenting cells; increasing responsiveness of cytotoxic lymphocytes to interferon-a; and decreasing frequency of regulatory and invariant T cells, including NKT cells and MAIT cells. Knowledge gained from the analysis of the neonatal immune landscape across GA provides a mechanistic framework to understand the unique susceptibility of preterm infants to both hyper-inflammatory diseases and infections

An immune clock of human pregnancy

The maintenance of pregnancy relies on finely tuned immune adaptations. We demonstrate that these adaptations are precisely timed, reflecting an immune clock of pregnancy in women delivering at term. Using mass cytometry, the abundance and functional responses of all major immune cell subsets were quantified in serial blood samples collected throughout pregnancy. Cell signaling-based Elastic Net, a regularized regression method adapted from the elastic net algorithm, was developed to infer and prospectively validate a predictive model of interrelated immune events that accurately captures the chronology of pregnancy. Model components highlighted existing knowledge and revealed previously unreported biology, including a critical role for the interleukin-2-dependent STAT5ab signaling pathway in modulating T cell function during pregnancy. These findings unravel the precise timing of immunological events occurring during a term pregnancy and provide the analytical framework to identify immunological deviations implicated in pregnancy-related pathologies

Multiomics modeling of the immunome, transcriptome, microbiome, proteome and metabolome adaptations during human pregnancy

Motivation Multiple biological clocks govern a healthy pregnancy. These biological mechanisms produce immunologic, metabolomic, proteomic, genomic and microbiomic adaptations during the course of pregnancy. Modeling the chronology of these adaptations during full-term pregnancy provides the frameworks for future studies examining deviations implicated in pregnancy-related pathologies including preterm birth and preeclampsia. Results We performed a multiomics analysis of 51 samples from 17 pregnant women, delivering at term. The datasets included measurements from the immunome, transcriptome, microbiome, proteome and metabolome of samples obtained simultaneously from the same patients. Multivariate predictive modeling using the Elastic Net (EN) algorithm was used to measure the ability of each dataset to predict gestational age. Using stacked generalization, these datasets were combined into a single model. This model not only significantly increased predictive power by combining all datasets, but also revealed novel interactions between different biological modalities. Future work includes expansion of the cohort to preterm-enriched populations and in vivo analysis of immune-modulating interventions based on the mechanisms identified. Availability and implementation Datasets and scripts for reproduction of results are available through: Https://nalab.stanford.edu/multiomics-pregnancy/

Integrated trajectories of the maternal metabolome, proteome, and immunome predict labor onset

Estimating the time of delivery is of high clinical importance because pre- and postterm deviations are associated with complications for the mother and her offspring. However, current estimations are inaccurate. As pregnancy progresses toward labor, major transitions occur in fetomaternal immune, metabolic, and endocrine systems that culminate in birth. The comprehensive characterization of maternal biology that precedes labor is key to understanding these physiological transitions and identifying predictive biomarkers of delivery. Here, a longitudinal study was conducted in 63 women who went into labor spontaneously. More than 7000 plasma analytes and peripheral immune cell responses were analyzed using untargeted mass spectrometry, aptamer-based proteomic technology, and single-cell mass cytometry in serial blood samples collected during the last 100 days of pregnancy. The high-dimensional dataset was integrated into a multiomic model that predicted the time to spontaneous labor [R = 0.85, 95% confidence interval (CI) [0.79 to 0.89], P = 1.2 × 10−40, N = 53, training set; R = 0.81, 95% CI [0.61 to 0.91], P = 3.9 × 10−7, N = 10, independent test set]. Coordinated alterations in maternal metabolome, proteome, and immunome marked a molecular shift from pregnancy maintenance to prelabor biology 2 to 4 weeks before delivery. A surge in steroid hormone metabolites and interleukin-1 receptor type 4 that preceded labor coincided with a switch from immune activation to regulation of inflammatory responses. Our study lays the groundwork for developing blood-based methods for predicting the day of labor, anchored in mechanisms shared in preterm and term pregnancies

Integrated plasma proteomic and single-cell immune signaling network signatures demarcate mild, moderate, and severe COVID-19

The biological determinants of the wide spectrum of COVID-19 clinical manifestations are not fully understood. Here, over 1400 plasma proteins and 2600 single-cell immune features comprising cell phenotype, basal signaling activity, and signaling responses to inflammatory ligands were assessed in peripheral blood from patients with mild, moderate, and severe COVID-19, at the time of diagnosis. Using an integrated computational approach to analyze the combined plasma and single-cell proteomic data, we identified and independently validated a multivariate model classifying COVID-19 severity (multi-class AUCtraining = 0.799, p-value = 4.2e-6; multi-class AUCvalidation = 0.773, p-value = 7.7e-6). Features of this high-dimensional model recapitulated recent COVID-19 related observations of immune perturbations, and revealed novel biological signatures of severity, including the mobilization of elements of the renin-angiotensin system and primary hemostasis, as well as dysregulation of JAK/STAT, MAPK/mTOR, and NF-κB immune signaling networks. These results provide a set of early determinants of COVID-19 severity that may point to therapeutic targets for the prevention of COVID-19 progression

Multiomics Characterization of Preterm Birth in Low- and Middle-Income Countries.

Importance: Worldwide, preterm birth (PTB) is the single largest cause of deaths in the perinatal and neonatal period and is associated with increased morbidity in young children. The cause of PTB is multifactorial, and the development of generalizable biological models may enable early detection and guide therapeutic studies. Objective: To investigate the ability of transcriptomics and proteomics profiling of plasma and metabolomics analysis of urine to identify early biological measurements associated with PTB. Design, Setting, and Participants: This diagnostic/prognostic study analyzed plasma and urine samples collected from May 2014 to June 2017 from pregnant women in 5 biorepository cohorts in low- and middle-income countries (LMICs; ie, Matlab, Bangladesh; Lusaka, Zambia; Sylhet, Bangladesh; Karachi, Pakistan; and Pemba, Tanzania). These cohorts were established to study maternal and fetal outcomes and were supported by the Alliance for Maternal and Newborn Health Improvement and the Global Alliance to Prevent Prematurity and Stillbirth biorepositories. Data were analyzed from December 2018 to July 2019. Exposures: Blood and urine specimens that were collected early during pregnancy (median sampling time of 13.6 weeks of gestation, according to ultrasonography) were processed, stored, and shipped to the laboratories under uniform protocols. Plasma samples were assayed for targeted measurement of proteins and untargeted cell-free ribonucleic acid profiling; urine samples were assayed for metabolites. Main Outcomes and Measures: The PTB phenotype was defined as the delivery of a live infant before completing 37 weeks of gestation. Results: Of the 81 pregnant women included in this study, 39 had PTBs (48.1%) and 42 had term pregnancies (51.9%) (mean [SD] age of 24.8 [5.3] years). Univariate analysis demonstrated functional biological differences across the 5 cohorts. A cohort-adjusted machine learning algorithm was applied to each biological data set, and then a higher-level machine learning modeling combined the results into a final integrative model. The integrated model was more accurate, with an area under the receiver operating characteristic curve (AUROC) of 0.83 (95% CI, 0.72-0.91) compared with the models derived for each independent biological modality (transcriptomics AUROC, 0.73 [95% CI, 0.61-0.83]; metabolomics AUROC, 0.59 [95% CI, 0.47-0.72]; and proteomics AUROC, 0.75 [95% CI, 0.64-0.85]). Primary features associated with PTB included an inflammatory module as well as a metabolomic module measured in urine associated with the glutamine and glutamate metabolism and valine, leucine, and isoleucine biosynthesis pathways. Conclusions and Relevance: This study found that, in LMICs and high PTB settings, major biological adaptations during term pregnancy follow a generalizable model and the predictive accuracy for PTB was augmented by combining various omics data sets, suggesting that PTB is a condition that manifests within multiple biological systems. These data sets, with machine learning partnerships, may be a key step in developing valuable predictive tests and intervention candidates for preventing PTB

Trip12, a HECT domain E3 ubiquitin ligase, targets Sox6 for proteasomal degradation and affects fiber type-specific gene expression in muscle cells

Abstract Background A sophisticated level of coordinated gene expression is necessary for skeletal muscle fibers to obtain their unique functional identities. We have previously shown that the transcription factor Sox6 plays an essential role in coordinating muscle fiber type differentiation by acting as a transcriptional suppressor of slow fiber-specific genes. Currently, mechanisms regulating the activity of Sox6 in skeletal muscle and how these mechanisms affect the fiber phenotype remain unknown. Methods Yeast two-hybrid screening was used to identify binding partners of Sox6 in muscle. Small interfering RNA (siRNA)-mediated knockdown of one of the Sox6 binding proteins, Trip12, was used to determine its effect on Sox6 activity in C2C12 myotubes using quantitative analysis of fiber type-specific gene expression. Results We found that the E3 ligase Trip12, a HECT domain E3 ubiquitin ligase, recognizes and polyubiquitinates Sox6. Inhibiting Trip12 or the 26S proteasome activity resulted in an increase in Sox6 protein levels in C2C12 myotubes. This control of Sox6 activity in muscle cells via Trip12 ubiquitination has significant phenotypic outcomes. Knockdown of Trip12 in C2C12 myotubes led to upregulation of Sox6 protein levels and concurrently to a decrease in slow fiber-specific Myh7 expression coupled with an increased expression in fast fiber-specific Myh4. Therefore, regulation of Sox6 cellular levels by the ubiquitin-proteasome system can induce identity-changing alterations in the expression of fiber type-specific genes in muscle cells. Conclusions Based on our data, we propose that in skeletal muscle, E3 ligases have a significant role in regulating fiber type-specific gene expression, expanding their importance in muscle beyond their well-established role in atrophy